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Abstract 

This is a sequel to the paper [K. Fujii : SIGMA 7 (2011), 022, 12 pages]. In this 
paper we treat a non-Gaussian integral based on a quartic polynomial and make a 
mathematical experiment by use of MATHEMATICA whether the integral is written 
in terms of its discriminant or not. 

1 Introduction 

The Gaussian is an abbreviation of all subjects related to the Gauss function ^-ip^^+i^+r) 
like the Gaussian beam, Gaussian process, Gaussian noise, etc. It plays a fundamental role 
in Mathematics, Statistics, Physics and related disciplines. 
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It is generally conceived that any attempts to generalize the Gaussian results would meet 
formidable difficulties. Hoping to overcome this high wall of difficulties of going beyond the 
Gaussian in the near future, a first step was introduced in [1], [2]. 

This is simply one modest step to go beyond the Gaussian but it already reveals many 
obstacles related with the big challenge of going further beyond the Gaussian. 

2 Cubic Case 

In this section we treat the cubic case. In the paper [ij^ the following "formula" is reported 



1) 



-D 

where D is the discriminant of the cubic equation 

ax^ + bx'^ + cx + d = 0, (2) 

and it is given by 

D = b^c^ + ISabcd - 4ac^ - 4b^d - TJa^d^. (3) 

The formula ([T]) is of course non-Gaussian. However, if we consider it in the framework of 
the real category then ([T]) is not correct because the left hand side diverges. In this paper 
we treat only the real category, and so a, 6, c, d, x, y are real numbers. 

Formally, by performing the change of variable x = tp, y = p for ([1]) we have 



LHS of© = j j e~'' \p\dtdp 



\a\e da I — = — =dt 
I ^(at3 + bt^ + ct + d)\ ^{at^ + bt^ + ct + d) 



by the change of variable a = \/ at^ + bt"^ + ct + d p. 
The divergence comes from 

J |cr|e~'^ da, 

^it is not easy for non-experts to understand this paper correctly 



while the main part is 

1 

-.dx 



^ [ax'^ + hx"^ + cx -\- d)\ ^ {ax^ + hx"^ + cx + d) 
under the change t — )■ x. As a kind of renormahzation the integral may be defined like 

X [ [ e-('^"'+''"'^+'="^'+'^^')(ia;(i?/ t = [ , ^ ^ dx. 

J Jr2 Jr I ^ (ax^ + bx"^ + cx + d)\^ {ax^ + bx'^ + cx + d) 

However, the right hand side lacks proper symmetry. If we set 

F{a,b,c,d)= [ [ e-K+'"'^+^"*''+''^')rfa;(i?/ (4) 
J Jdr 

where Dr = [—R, R] x [—R, R] (i? ^ 0), then it is easy to see 

F{—a, —b, — c, —d) = F{a, b, c, d). 

Namely, F is invariant under Z2-action. This symmetry is important and must be kept even 
in the renormahzation process. The right hand side in the "definition" above is clearly not 
invariant. Therefore, by modifying it slightly we reach the renormalized integral 

Definition 

t [ [ e-(''"'+'"'^+^"J^'+'^^')(ixrf2/ t = [ , ^ dx. (5) 

J JR2 Jr \/ {ax^ + bx"^ + cx + dy 

The Gamma-function T{p) is defined by 

POO 

T{p) = / e'^'xP^^dx {p > 0) (6) 
Jo 

and the Beta-function B{p, q) is 

B{p,q)= [ xP-\l-xy-^dx {p,q>0). (7) 
Jo 

Then the result in [2] is 

Formula 
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(I) For < 



where 



C_ = v^fi(-,-). 



(II) For D > 



where 



R (ax^ + + cx + rf)2 v^D 



(Ill) C_ and C+ are related by C+ = -\/3C_ through the identity 

v3B{y-)=nB{y-). (10) 

See [2] in detail. Our result shows that the integral depends on the sign of D. This 
formula has been conjectured by Morozov and Shakirov [1] in a different context. 

A comment is in order. If we treat the Gaussian case (: Q-i"-^^+b^y+<^y^)"^ then the integral 
is reduced to 



i 



1 ,,^j^^md> (111 



In ax2 + bx + c \J-D V-D 
if a > and D = — Aac < 0. Because 

^ r(i)r(i) 1 1 

1 r(i) ^2' 2^ 

Let us check whether the renormalized integral ^ is reasonable or not by making use of 
the results. 

For the integral (j4]) it is easy to see 

^^-^^)i^(a,6,c,d) = j 1^ (x3./_^2^.^y2)^-(-3+^-^.WW)d^rfy = 0, 

^^-^^)F(a,6,c,d) = f (x2y.x2y-x3.xy2)e-(-W.WW)d^dy = 0, 

(12) 

(|^-^^)na,6,C,d) = j 1^ (xy2.^y2_^2^.^3)^-(a.3+,..,WW)rf^rfy = 0. 
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On the other hand, if we set 

J^ia, b, c,d)= [ — dx = — 

Jr ^(ax3 + 6x2 + CX + df V±D 

then we can also verify the same relations 

[d-add-dbd-c) ^Mc,d) = 0, 
f d d d d\ 

f d d d d\ 

[d-cd-c-mdd) -^(-'^'^'^) = o- 

by use of MATHEMATICA. 

Therefore we can say that the definition ([5]) is not so bad (maybe, good). 



3 Quartic Case : Mathematical Experiment 

In this section we treat the quartic case. The discriminant of the quartic equation 

ax^ + bx^ + cx"^ + dx + e = (14) 

is given by 

D = 256a^e^ - 4b^d^ - 27a^d^ - 276^e^ - 1280^0^6^ + b'^c^d^ + IGac^e 
-Aac^d^ - ib^c^e + 1440^0^^6 - 6ab^d^e + lUab^ce^ - 192a^bde^ 
+18abcd^ + ISbhde - SOabc^de. (15) 

See for example [2]. 

Here we consider a non-Gaussian integral 

= T{a, b, c,d,e) = J J e-{--'+bx'y+cx^y'+d^y'+ey^) ^xdy (16) 

and study whether this integral can be written in terms of its discriminant (11 5p or not like the 
formula in the preceding section. The conclusion is negative, while we have very interesting 
mathematical "phenomena" stated in the following. 
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By use of the same change of variable x = tp, y = p in the cubic case we have 
:F= j j e~^'(^''+'"'+'^''+''+')\p\dtdp = j y e~(^''+''''+'''"+''+'y\p\dp^dt 
and the change of variable a = \J at^ + bt^ + ct"^ + dt + ep {at^ + bt^ + ct'^ + dt + e > 0) gives 



^4 



J- = / lale^^'^da / ^ ^dt 

^{at^ + bt^ + ct2 + dt + e)2 

/tt /■ 1 

=dx (17) 



2 Jr -y/ (ax"^ + + cx^ + (ix + e)^ 



because 



OO /"OO 



Problem I How can we calculate f[T7|) 



? 



There is no method to calculate at the present time, so we make a mathematical ex- 
periment by use of MATHEMATICA. From the lesson of quadratic and cubic cases it may 
conjecture 

^=^^ (D<0, a>0) (18) 
where C is some constant and 12 = 4 x 3. 

For the integral of exponent (|T6|) we have a system of differential equations 

= 0, 

= 0, 

J" = 0, 

= 0, (19) 

= 0, 

= 0. 
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The proof is straightforward. For example 
f d d d' 



da dc db' 
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On the other hand, if we set 



then we have 



f d d 
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~ _ (c2-3M+12ae)(3(i2-8ce) 
J~ — — (_/ 



J- = C 



f = C 



36(-D)i3/i2 

^c^ -3M+ 12ae)(c(i-66e) 
18(-D)i3/i2 ' 

'c' - 3bd + 12ae){c' - 2bd - 4ae) 



9(-D)i3/i2 
c' -3bd+12ae){16ae-bd) 



F = -C- 



36(-D)i3/i2 
3M+ 12ae)(6arf- 6c) 



18(-D)i3/i2 
c' - ?,bd + I2ae){?,b' - ^ac) 



36(-D)Wi2 

by use of MATHEMATICA (verification by hand is very tough). 

As a result T T . However, from (^11) we obtain an interesting quantity 

E = c' — 3bd + 12ae ax^ + bx^ + cx' + dx + e. 

It is not clear at the present time what £" is, so we present 
Problem II Make the property of E clear. 
We believe that E will play an important role in the calculation. 
In last, we note some interesting fact. If we set 

(O>0) 



then the same relations fl2Tl) hold apart from the sign 
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^ _ Jc'^ -3bd+12ae){3d'^ -See 



T = -C- 



c'^ -3bd+12ae){cd-6be) 



T = C 



c 



132)13/12 

3bd+12ae){c^ - 2bd - Aae) 
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7 = -C 
T = C 



{c^ -3bd+12ae){16ae-bd) 
36D13/12 

c2 - 3bd + 12ae){6ad - be) 



7 = C 



18^13/12 

(c2-3M+12ae)(3b2-8ac) 



36DW12 

These "phenomena" are interesting enough and worth studying in detaiL 



4 Concluding Remarks 

In this note we treated a non-Gaussian integral based on a quartic polynomial and made a 
mathematical experiment by use of MATHEMATICA. Though our work is far from obtaining 
an explicit value of the integral, we found an interesting quantity E. In order to make it 
clear hard work will be needed. 

See [3], |1] as recent results on this topic and as a general introduction to non-linear 
algebras. We expect young mathematicians or mathematical physicists to take part in this 
fascinating topic. 
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